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ACOUSTIC  RADIATION  AND  SCATTERING  FROM  ELASTIC  STRUCTURES 

PRKIS 

1 • This  report  presents  a.  numerical  technique  for  the  linear  dynamic  analysis 
^ finite  elastic  structure  immersed  in  an  infinite  homogeneous  acoustic  medium* 
It  is  required  to  determine  the  vibratory  motion  of  the  structure  and  also  the 
associated  acoustic  field  in  the  fluid,  when  the  structure  is  either  subjected 
to  internal  applied  forces  or  is  acting  as  a scatterer  of  an  incident  acoustic 
wave,  A finite  element  analysis  of  the  structure  is  matched  at  the  structure- 
fluid  interface  with  an  integral  equation  representation  of  the  exterior  acoustic 
field,  leading  to  a coupled  system  of  equations  which  may  be  cast  in  either 
acoustic  or  structural  form*  The  former  approach  is  preferred  here  for  which 
numerical  results  are  presented  when  the  method  is  applied  to  plane  wave  scattei'- 
ing  by  thick  and  thin  elastic  spherical  shells* 

CONCLUSIONS 

2*  By  combining  equations  derived  from  a finite  element  analysis  of  a vibrat- 
ing elastic  structure  with  those  from  an  integral  equation  representation  of  an 
infinite  exterior  acoustic  field,  an  analysis  of  the  coupled  dynamic  interaction 
problem  has  been  shown  to  be  feasible* 

3*  The  structural  and  acoustic  models  may  be  chosen  completely  independently 
although  the  procedure  does  simplify  if  inteipolation  nodes  are  made  to  coincide, 
the  coupling  being  defined  through  three  interaction  matrices  X,  L and  L»*  Thus 
advantage  may  be  taken  of  existing  computer  program  packages  designed  to  solve 

the  two  uncoupled  problems*  In  particular,  it  is  hoped  to  combine  the  PAPEC^^ 
structural  analysis  program  with  the  acoustic  radiation/scattering  program  des- 
cribed in  this  report^®* 

4*  Some  test  examples  for  which  analytic  solutions  are  available  have  demon- 
strated the  practicability  and  accuracy  of  the  approach*  It  is  hoped  in  the  near 
iuture  to  be  able  to  con5)are  numerical  and  experimental  results  for  some  more 
complex  structiires  in  water*  If  an  unacceptably  large  number  of  degrees  of  free- 
dom for  the  structural  model  are  required,  it  may  then  be  necessary  to  consider 
the  alternative  modal  approach* 


5.  This  report  presents  a numerical  technique  for  the  linear  dynamic  analysis 
of  a finite  elastic  structure  immersed  in  an  infinite  homogeneous  acoustic  medium. 
This  problem  commonly  occurs  in  underwater  acoustics  where  it  is  of  considerable 
interest  to  determine  the  acoustic  field  both  radiated  by  a submerged  vibrating 
structure!  and  also  scattered  by  a submerged  elastic  structure.  It  is  well-known 
that  the  vibrational  properties  of  a structure  can  significantly  effect  the 
scattered  acoustic  field|  particularly  when  the  acoustic  medium  is  wateri  the 
impedance  mismatch  being  much  less  than  between  a structure  and  air,  and  the 
assumption  which  is  frequently  made  that  the  structure  is  perfectly  rigid  is  often 
an  oversimplification  of  the  true  situation. 

6.  Analytical  approaches  to  such  coupled  structure-fluid  interaction  problems 
eure  almost  invariably  concerned  with  spherical  or  infinite  cylindrical  geometries 
for  which  the  classical  method  of  separation  of  variables  is  available^  since 

in  these  cases  normal  structural  modes,  not  coupled  by  radiation  loading,  exist. 

In  particular,  approximate  shell  theory  has  been  extensively  used  to  determine 

the  response  of  thin  elastic  shells  under  acoustic  loading.^  ^ 

7.  Until  relatively  recently  little  progress  had  been  made  towards  the  solu- 
tion of  interaction  problems  for  other  geometries.  However,  now  that  general 
numerical  methods  have  been  independently  developed  for  both  the  structural  and 
acoustic  problems,  the  feasibility  of  combining  these  techniques  to  tackle  the 
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coupled  problem  has  been  realised.' 


For  a complex  stiucture  subjected  to  known 
.13 


applied  forces  the  finite  element  method  has  become  an  accepted,  well-proven, 
and  highly  successful  analysis  tool.  Similarly,  although  certainly  less  widely 
applied,  methods  for  determining  the  acoustic  field  radiated  or  scattered  from  a 
structure  have  been  developed  when  boundary  conditions  on  the  structure  surface 
are  assumed  to  be  known,  e.g.  rigid  surface,  soft  'pressure  release'  surface. 

These  methods  are  usually  based  on  integral  equation  formulations^ and 
^proximated  using  finite  element  type  expansions. 

8.  Indeed,  although  the  term  'finite  elements'  is  only  now  being  explicitly 
associated  with  methods  for  the  solution  of  integral  equations,  if  the  name  is 
understood  to  refer  to  the  local  piecewise  nature  of  the  approximation  to  the 
domain  of  the  equation  and  also  to  the  unknown  fxmction,  and  not  necessarily 
cissociated  with  any  particulsir  method  for  determining  the  function  parameters, 
such  techniques  have  been  in  evidence  for  some  time.  Most  commonly,  when  apply- 
ing the  finite  element  method  to  integral  equations,  the  parameters  are  deter- 
mined through  collocation  since  the  more  familiar  methods  of  Galerkin  or  using  a 
variational  principle  require  a further  integration  over  the  domain  of  the 
equation.  When  that  domain  is  an  arbitrary  surface  in  three-dimensional 
Euclidean  space  over  which  the  integrations  are  performed  numerically,  a repeated 
integration  could  be  extremely  costly  in  computer  time.  In  fact  Maxwell  effec- 
tively solved  the  integral  equation  relating  the  potential  of  a thin  charged 
conducting  sq[uare  plate  to  the  s\xrface  charge  density  by  a finite  element 
1 8 

e^proach  in  1879  • Although  the  complete  details  are  not  given  he  apparently 

accounted  for  the  expected  singularity  in  the  charge  density  at  the  edges  and 
corners  of  the  plate  by  using  singular  basis  (shape)  functions  in  these  regions. 

19 

Also  product  integration,  which  was  introduced  by  Young  in  1954  as  a technique 
for  the  solution  of  integral  equations,  is  in  its  most  popular  form  no  more  than 
a finite  element  method. 
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ACOUSTIC  EQUATIONS 


9«  SmaJ.1  amplitude  acoustic  wraves  propagate  through  an  ideal  homogeneous 
fluid  of  density  and  speed  of  sound  o,  according  to  the  linear  wave  equation 

P(r,t)  = 0 (1) 

c 3t 

where  P(r,t}  is  the  excess  acoustic  pressure  at  a position  r in  the  fluid  at 
time  t,  Sound  pressure  is  related  to  particle  velocity  V(r,t)  through 

3V 

VP(£,t)  = -Pf  ^ Cr,t).  (2) 

iwt 

ID,  For  a single  frequency  harmonic  time  dependence  of  the  form  e , where  u 
is  the  anguliir  frequency,  equation  (l  ) becomes  the  Helmholtz  or  reduced  wave 
equation 

(V^  -f  k^)  p(r)  = 0 (3) 

where 

P(r,t)  = p(r)e^"^  (if) 

and  k = w/c  is  the  acoustic  wavenumber. 

11.  When  the  domain  of  the  partial  differential  equation  (3)  is  that  infinite 
region  exterior  to  the  structure  surface,  techniques  normally  employed  for  the 
approximate  solution  of  such  equations  (finite  differences,  finite  elements  etc) 
present  obvious  difficulties  of  implementation,  A boundary  must  bo  introduced 
at  sc*ne  finite  distance  from  the  structure  together  with  a boundary  condition 
there  which  should  approximate  the  true  radiation  condition  at  infinity,  to  ensure 
that  al]  aco’jstic  waves  either  radiated  or  scattered  by  the  structure  are  outgoing 

20 

there.  Such  an  approach  was  suggested  by  Zienkiewicz  and  Newton  using  a system 
of  dashpots  on  the  outer  fluid  boundary,  although  the  correctness  of  this  pro- 

21 

cedure  has  been  questioned.  Even  if  the  appropriate  boundary  condition  could 
22 

be  derived  , such  methods  necessarily  require  a large  number  of  nodes,  particu- 
larly in  a three-dimensional  situation,  at  each  of  which  an  approximation  to  the 
solution  is  produced.  In  practice  however,  it  is  often  only  the  acoustic  field 
at  a selected  number  of  nearl’ield  positions  and  also  the  farfield  radiation 
pattern  that  is  of  interest. 

12,  Thus  the  majority  ol'  workers  in  acoustics  have  chosen  to  reformulate  the 
dil'ferent ial  problem  (i.e.  differential  equation  plus  sufficient  boundary  condi- 
tions' a • in  integral  equation.  This  approach  has  the  immediate  advantage  that 
the  ini'j  nite  exterior  domain  of  the  differential  equation  reduces  to  a finite 
douwin  of  one  dimension  less  (the  structure  siirface)  for  the  integral  ec^ation. 
Hence  only  surface  values  are  initially  calculated  from  which  the  field  at  any 
posit..ons  ol  interest  may  be  evaluated  through  the  integral  representation 
inherent  l.i  the  approach.  Also,  the  exact  radiation  condition  is  automatically 
inciadbd  m sny  integral  equation  formulation. 


13.  Unfortunately,  all  the  classical  formulations  of  the  Helmholtz  equation 
in  an  exterior  jregion  as  an  integral  equation  ■■  equivalent  single  or  double 
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layer  source  distributions  cr  via  Green's  Theorem  - break  dovm  either  through 
non-'uniqueness  or  non-existence  of  their  solutions  at  certain  wavenumbers. 

This  failure  is  completely  non-physieal  in  nature,  being  solely  due  to  the 
integral  equation  formulation,  and  is  associated  with  the  existence  of  eigen- 
solutions  of  the  Helmholtz  equation  within  the  region  interior  to  the  surface 

of  interest.  * Various  'improved'  formulations  have  been  proposed  which  are 

i L.  9R 

either  designed  to  mitigate  the  problem  at  low  frequencies,  ^ or  to  eliminate 
26*28 

it  completely.  These  latter  methods  require  considerably  more  computa- 

tional effort  to  implement  and  have  as  yet  not  f oimed  the  basis  of  any  practical 
acoustic  radiation/scattering  computer  program,  although  test  programs  have  been 

29  30 

written  and  applied  to  some  simple  problems  in  two  and  three  dimensions.  A 
comprehensive  feasibility  study  for  such  a program  package  has  recently  been 
conqjleted."^^ 

14.  In  order  to  Illustrate  the  coupling  procedure  proposed  here,  the  integral 

14 

equation  approach  due  to  Schenck  is  employed,  edthough  within  the  possible  fre- 
quency limitations  previously  mentioned  any  available  integral  equation  formula- 
tion could  be  used. 

15.  Consider  an  acoustic  wave  (Pj^q(£))  incident  upon  the  closed  surface  S of 

the  elastic  structure,  resulting  in  a scattered  acoustic  wave  with  pressure 
P ...(n).  Then,  from  Green's  Theorem,  the  total  acoustic  pressure 

pCO  (®  ^scatt^^^^  satisfies  the  Helmholtz  Integral  formulae"^ 


■ Pino^£‘^ 


r'  in  E (5a) 


/ {p(£)  f-  (£.£■)  - - 'in=‘£’)  £'  »"  = (^i.) 

1 - Pinc^£'^  ^ 

where  B is  the  infinite  acoustic  medium  exterior  to  S,  1 is  the  interior  of  the 
surface  S and  a is  the  solid  angle  subtended  by  the  acoustic  fluid  at  the  point 
on  S (for  a anooth  surface,  a = 2v  everyidiere). 


G(r,r')  = 


-ik  r-r' 
e ^ ' 


is  the  three-dimensional  free-space  Green's  Ainction  for  the  Helmholtz  equation 
(5),  and  — p denotes  differentiation  in  the  outward  normal  direction  (n)  at 

X* 

^ on  S (i.e.  from  S into  the  fluid).  In  the  case  of  pure  radiation  the  total 

acoustic  field  p(r),  which  now  consists  of  radiated  pressure  only,  also  satisfies 

the  Helmholtz  formulae  (5)  with  the  term  p.  (r)  absent. 

^ino  — ' 

1 6.  On  any  part  of  the  surface  S viiich  has  a well-defined  normal  direction, 
the  relation  (2)  gives 

|j^(r)  = -Iwpj.  v(r)  (7) 

where  v(r)  is  the  outward  normal  surface  velocity. 
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17.  Hence,  equations  (;.>b)  and  (be),  known  as  the  Sui'faco  Helmholtz  Equation 
(she)  and  interior  Helmholtz  Relation  (IHR)  respectively,  relate  the  acoustic 
pressure  on  the  structure  surface  to  the  normal  surface  velocity.  Both  have 
been  used  as  the  basis  for  the  solution  of  acoustic  radiation  or  scattering 
probltvns  «*ien  either  the  surface  pressure  or  velocity  is  known,  although 

neither  is  entirely  adequate. 3h>''b  most  ajitable  for 

numerical  treatment  but  exMbits  non-uniqueness  of  the  solution,  whenevo'  the 
wavenumber  k is  an  eigenvalue  of  the  Dirichlet  eigenproblem  for  the  interior 

region  I The  i“R  can  be  shown  to  possess  a unique  solution  if  it  is  satis- 

35 

fied  at  all  interior  points'^  (for  an  axisymmetric  problem  it  is  only  necessary 
to  enforce  tiie  relation  at  all  points  along  the  axis  of  symmetry)  but  suffers 
from  numt-rical  instability  due  to  its  similarity  to  an  integral  equation  of  the 
first  kill!  with  a non-singular  kernel. 

1S.  The  idea  of  Schenck  was  to  combine  the  numerical  stability  of  the  SHE  with 
the  ur.ioueness  property  of  the  IHR  by  overdetermining  the  set  of  linear  equations 
obtained  from  a numerical  approximation  of  the  SHE  with  a small  number  of  addi- 
tional equations  derived  from  the  IHR.  The  resulting  linear  system  is  then 
solved  by  a least  squares  procedure.  This  method  is  generally  known  by  the 
acronym  CHIEF  (Combined  Helmholtz  I^ntegral  Equation  formulation).  Schenck  proved 
that  at  a critical  wavenumber  only  the  required  solution  of  the  SHE  will  simul- 
taneously satisfy  the  IHR,  provided  the  IHR  is  not  evaluated  on  a nodeil  surface 
of  the  relevant  interior  eigenfunction. 


19.  In  practice  the  CHIEF  method  is  satisfactory  only  at  relatively  low  fre- 
quencies when  a sufficient  number  of  correctly  positioned  interior  nodes  can  be 

2S 

selected  if  required.  Jones  ' has  suggested  a systematic  approach  to  this  choice. 

20.  In  order  to  obtain  a discrete  matrix  af)proximation  to  the  SHE  and  IHR  a 
fi.iite  element  method  is  used.  Approximations  to  the  surface  acoustic  pressure 
and  ncrmal  surface  velocity  are  assumed  in  the  forms. 


m 

p(r)  L 
i=1 


v(r)  Z 
i=1 


'i 


(8) 

(9) 


where  and  Af^Cf)  are  scalar  basis  Ainctions  defined  on  some  surface  S* 

approximating  S.  Collocation  of  the  SHE  (5b)  at  m surface  points  r . , 

iri  ,m,  at  which  there  exists  a well  defined  normal  direction,  and  the 

IHR(bc''  at  m'  interior  points  i=m^1  , m+2,...,  m>m'  leads  to  the  matrix 
system 


A?  = Bv  > p 


inc 


wher  0 , 


(10) 


(11) 
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for  1=1,2,...,  iiH-m';  J=1,2,...,  m, 

c T f T 

and  p = Ip^  P2  ...  Pjjjl  , V = {v^  Vg  ...  , and 

Pino  = •Pino<J:i>  Pino‘i2>  — Plno^ii^.'' 

Details  of  a particular  choice  of  basis  functions  and  a method  of  generating 
the  approximate  surface  S*  are  discussed  later. 

21.  Note  that  the  resulting  matrices  are  full,  in  contrast  with  the  sparse 
banded  matrices  that  are  typical  of  the  numerical  solution  of  differential 
equation  problems  by  finite  element  or  finite  difference  methods,  and  the  coef- 
ficients, which  are  frequency  dependent,  must  in  general  be  derived  via  time- 
oonsiimlng  numerical  quadrature.  Also,  since  collocation  has  been  employed  to 
obtain  sufficient  relations  between  the  parameters,  the  system  is  not  symmetidc. 

22.  The  system  of  equations  (l0)  is  sufficient  to  determine  completely  the 
acoustic  field  in  two  particular  oases  of  interest;  namely  for  radiation 
(Pj^^(jp)  = 0)  where  the  surface  velocity  is  specified,  and  for  scattering  irtiere 

the  acoustic  Inqpedance  of  the  surface  is  specified,  most  usually  assumed  to  be 
a perfectly  hard  or  rigid  surface  (v(r)  *0), 

23.  However,  in  general  the  vibrational  response  of  the  strvicture  must  be 
taken  into  account,  in  idiloh  case  equation  (IO)  merely  provides  a relationship 
between  acoustic  pressure  and  normal  velocity  at  the  structure  surface. 

STRUCTURAL  EQUATIONS 

24.  The  region  enclosed  by  the  surface  S is  assumed  to  contain  an  elastic 
structure  whose  motion  is  governed  by  the  linear  equations  of  elasticity, 

a2 

div  a(r,t)  + F(r,t)  = p — w U(r,t)  (l3) 

OX 

where  o{r,t)  is  the  stress  tensor,  £(r,t)  represents  external  forces,  is  the 
structure  density,  and  U(r,t)  is  the  particle  displacement. 

25.  The  discretisation  of  these  equations  by  the  finite  element  method  is  now 

a familiar  technique^ ^ and  will  not  bo  described  in  detail  here.  An  approxima- 
tion to  the  structural  displacement  is  assumed  in  the  form 


U(r,t)  = 


J,  Jii  ‘i(£) 


where  each  is  a diagonal  matrix  (of  order  d,  the  spatial  dimension  of  the 

problem)  of  basis  functions  defined  throughout  the  volume  of  the  structure. 
These  basis  functions  are  normally  chosen  such  that  they  vanish  at  all  but  one 
node  of  an  element  subdivision  of  the  structure,  in  iriilch  case  the  vector 


a. 

parameter  describes  the  displacement  at  that  node*  If  each  component  of  the 
parameter  has  the  same  associated  basis  function,  as  Is  frequently  the  case 
In  priu'luf,  the  matrix  may  be  replaced  by  a scalar  basis  function  ^^(r). 

2fa»  The  finite  element  equations  for  the  structure  with  the  harmonic  time 
dependence  omitted,  are  then  of  the  form 

(K  iwC  - w^M)  q = f^'^^  ♦ f^P^  (15) 

where  K,  C and  M are  stifl’noss,  damping  and  mass  matrices  respectively,  q is  the 

(k) 

vector  of  parameters  u^,  f'  is  a consistent  load  vector  derived  from  known 

applied  forces,  and  f(p)  la  u consistent  load  vector  representing  the  acoustic 
fluid  pressure  acting  on  the  fluid-atruoturo  boundary.  Explicitly,  the  vector 

f^^^  has  components 

where  this  notation  is  used  to  define  simultaneously  all  the  components  of  f 
associated  with  the  vector  parameter  \^* 


2/,  The  matrices  K and  M and  usually  C are  synmetrio  and  baitded  and  if  struc- 
tural danqaing  is  ignored,  the  system  is  re«a.  If  there  is  no  acoustic  loading 
on  the  structure  or  if  the  structure  is  only  subjected  to  a static  fluid  pressure, 
the  sy’.+cj;  jf  equations  (l5)  alone  defines  the  displacement  of  the  structure  in 
terms  of  the  applied  loads.  However,  in  the  case  of  dynamic  fluid-solid  inter- 
action, these  equations  only  relate  structural  displacements  to  the  aooustic 
px-essure  at  the  fluid-struotui'e  interface. 

COUPLED  SgUATIONS  OF  MOTION 

26.  The  complete  solution  of  the  fluid-structure  Interaction  problem  may  now  be 
described  by  combining  the  acoustic  equations  (l0)  with  the  stiructural  equations 
(15).  Substitution  of  the  representation  (8)  for  the  surface  aooustic  pressure 
into  (16)  leads  to 


= -Lp 


where 


(17) 


0j(r)  *.(r) 


n dC 

r 


(18) 


for  i=1,2,  ...,  n;  J=1  ,2,  ...,  m,  which  ensures  continuity  of  sound  pressure 
between  the  structural  and  acoustic  models. 


29.  Continuity  of  normal  surface  velocity  is  achieved  by  matching  at  some  set 
of  txurface  nodes,  r|,  i=1  ,2,  ...,  k,  the  approximations  (9)  and  (14)  to  give 

Xv  iwL'q 


(19) 


wher« 


‘ij  ‘ (2®) 

^ OX*  •••f  ID|  And 

■-'ij  ■ ‘j‘2:i)i!l  (»') 

for  lei ,2,  ...,  k;  Jel ,2,  ...,  n,  with  n*  baln^  the  outward  nortaal  direotion 

at  r!  • 

— i 

30,  Normally  the  baala  functions  <p.{r)  and  y^Cr)  would  be  chosen  to  bo  the 

same  and  to  Interpolate  m surface  points.  These  would  then  be  the  natural 
collocation  points  for  the  SHE  and  would  also  be  the  positions  at  which  to 
ensure  the  continuity  of  velocity.  In  this  case  X would  be  the  identity  matrix. 
If  these  acoustic  nodes  coincide  with  some  or  all  of  the  surface  structural 
nodes  then  the  non-sero  elements  of  the  matrix  L*  are  simply  components  of  the 
outward  normal  dlreotlons  at  the  acoustic  nodes. 

31.  The  solution  of  equations  (10),  (15),  (l 7)  and  (l9)  may  be  accomplished  In 
a number  of  ways.  Elimination  of  the  structural  displacement  vector  leads  to  a 
combined  matrix  equation  In  acoustic  form 

(A  ♦ DL)  p = Df^*'^  ♦ p,_  (22) 


where 


D = iwBX'^L'  (K  ♦ i«C-«Sl)"\ 


Once  the  surface  (U'essure  has  been  determined  the  structural  motion  follows 
directly  from 


q » (K  ♦ 1«jC  - wSl)"'*  (f^''^  - Lp). 


Initial  elimination  of  the  acoustic  pressure  leads  to  a perturbation  of  the 
structural  equations, 

(K  ♦ iwC  - wSl  ♦ IwU'^BX'^L')  q = f^*'^  - U“^  p^^^  (25) 

witli  the  surface  pressure  then  being  toy 

p = A ^ (u'BX”'l'  q ♦ Pjjjp)- 

yi.  If  it  is  necessary  to  I'ona  non-squ<ire  m a trices  A and  B to  overcome  non- 
uniqueness problem:',  the  inverse  of  matrix  A Is  to  be  understood  In  a i^eneralised 
least  squares  sx'rux'. 


U’'bx“’l'  = ^(R^^IK^) 


then  equation  (25)  becomes 


Jr 


ri ' 


the  fluid  la  oonsl dared  Inocmpresalble  (the  low  frequency  limit,  k c o)  and 
hence  unable  to  auatain  a aound  field,  the  matricea  A and  B become  totally  real 
and  Ima^nary  reapeotively  and  the  added  damping  term  ia,  aa  expected,  iden- 
tically aero.  In  thla  oaae  it  la  the  Laplace  equation  which  governa  the  fluid 
behaTiour,  for  which  there  are  no  non-unlqueneaa  or  non-exiatenoe  problema 
aaaoolated  with  an  integral  equation  Asrmulation. 

33*  Krom  a oomputatlonal  point  of  view  the  aystem  of  equationa  (22)  la  denae, 
complex,  unaymmetrlo  and  of  dimenaion  (m^m*)  x m.  Aa  an  Intermediate  atep  it 

(k ) 

la  neceaaary  to  form  the  producta  DL  amd  Df'  , not  by  computing  explicitly  the 

Inver ae  matrix  (K  * iuC  - but  by  aolvii^  an  nd  x nd  ayatem  of  linear 

equationa  with  (K  ♦ iwC  - wSi)  aa  the  coefficient  matrix  and  the  columna  of  L 
(k) 

and  alao  f'  * aa  the  right  hand  aide  veotora.  Efficient  routinea  are  available, 
due  to  the  aymmetry  and  bandedneaa  of  the  coefficient  matrix,  to  carry  out  thia 
procedure  the  reault  of  which  ia  retained  for  aubaequent  evaluation  of  the 
atruotural  dlaplaoement  parametera. 

y*..  Now  conaider  the  ^atem  of  equationa  (23}«  Due  to  the  aoouatio  perturbing 
term,  the  ayatem  no  longer  poaaeaaea  the  atruotural  featurea  of  bandedneaa  and 
ayrmaetry,  but  ia  relatively  denae,  oonplex,  unaymmetrlo  and  of  dimenaion  nd  x nd. 

Aa  an  Intermediate  atep  in  thla  oaloulatlon  it  ia  neceaaary  to  form  A and 

A solving  an  (m  ♦ m* ) x m denae  complex  ayatem  with  multiple  right  hand 

sidea  uaing  a leaat  aquarea  procedure. 

33.  Since  in  general  n » (m  -f  m')  it  would  appear  computationally  more  effi- 
cient to  formulate  the  coupled  aolution  aa  the  ayatem  of  equationa  (22)  and  thla 
la  preferred  here,  although  aome  authora  have  adopted  the  alternative 

. 8,12 

approach. 

3b.  For  many  real  complex  atruoturea  aubmerged  in  water,  particularly  when  no 
synmetry  ia  preaent,  the  finite  element  model  would  neoeaaarily  have  a large 
number  of  degreea  of  freedom  and  thua  exoeaaive  computer  coat  may  reatrlot 
eithor  of  the  above  approaohea.  In  thia  oaae  a modal  approach  may  be  uaed  to 

advantage, whereby  the  atructure  dlaplaoement  ia  ^proxlmated  by  a linear 
combination  of  the  dominant  In-vaouo,  undamped  normal  modea  of  the  atructure. 

The  modal  frequenciea  (a^)  and  mode  ahapea  {6^)  are  efficiently  obtained  from 

equation  (l3)  with  the  load  and  damping  terma  aet  to  aero. 


(K 


A^  M)  = 0. 


A convenient  normal iaation  ia  uaually 
1 


(29) 


(30) 


T 

®i  «1 


(31) 


T 


1 

I 


■ 1 
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with  the  orthof^ona.!! ty  relations 

. 0 (32) 

holding  for  i ^ J. 

If  the  atructvire  dis  pla oeaent  is  now  expressed  in  terms  of  N of  these  modes  as 

q = 0 a (35) 

where  a is  a vector  of  parameters  and  6 is  the  natrix  'xrhose  oolumns  are  the 
N normal  modes,  the  aooustlo  and  structural  forms  of  the  coupled  problem  become 
respectively, 

(A  ♦ SL)  p = Sf^*"^  + p^^^  (34) 

where 

S = iwBX’^L'  9 (a  + iw  e^C  6 - w^l)"''  9^  (35) 

and 

(A  ♦ iu.  (C  ♦ u"‘'bL')  9 - Jl)  a . 9'^(f^‘'^  - 1^"%,  ),  (36) 

ino 

2 

A is  the  diagonal  matrix  with  elements  • 

37.  These  systems  are  similar  to  those  given  previously  except  that  if  only  a 
small  number  of  dominant  elgenmodes  are  chosen  to  represent  the  displacement, 
the  matrices  derived  from  the  structural  model  are  of  smaller  size  (N  « n).  Of 
coui’se  the  eigensolutions  of  the  structure  also  have  to  be  determined,  but  this 
need  only  be  caiTied  out  once  for  a structure  while  a coupled  solution  might  be 
required  at  a number  of  frequencies.  If  no  internal  damping  is  present  or  if 
present  can  be  represented  in  a particularly  simple  manner,  the  matrix  inversion 
required  in  equation  (35)  becomes  trivial, 

56.  Once  the  surface  pressure  and  structural  displacement  (and  hence  normal 
siufaoe  velocity)  have  been  evaluated,  the  pressure  at  any  position  in  the  fluid 
may  be  determined  directly  from  the  Helmholtz  equation  (5a)  via  numerloal 
quadrature  in  an  analogous  manner  to  the  derivation  of  the  coefficients  of  the 
acoustic  matrices  A and  D.  For  positions  at  large  distances,  R,  in  terms  of  the 
acoustic  wavelength,  from  the  structure  the  radiated  or  scattered  part  of  the 
acoustic  field  tuis  an  angular  distribution  independent  of  R.  The  dependence  on 
K is  of  the  form  fixp(-lkR)/Jl  irtiloh  may  be  factored  out  of  the  Helmholtz  integral 

58 

to  simplify  its  evaluation. 

y.K  Thus  equations  which  describe  a complete  solution  of  the  coupled  fluid- 
structure  interaction  problem  have  been  presented  for  a structure  either  subjected 
to  internal  applied  forces  or  acting  as  a soatterer  of  an  incident  acoustic  wave. 

(k) 

Noimaily  one  of  the  vectors  f'  ' and  p^^^^  will  be  identically  zero  although  the 
formulation  does  not  require  this  to  bo  so. 


f 
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NUMHIICAL  IlCPLEHBWrAriOW  AND  RBSUIiTS 

ifO*  In  or^eap  to  evaluate  the  ooeffiolenta  of  the  satrloes  A and  B in  the 
aoouatlo  equations  (10)  it  is  neeessaiy  to  define  a surface  S*,  an  approximation 
to  S,  choose  the  basis  functions  and  x^ir)  appearing  in  the  surface  pres- 

sure and  normal  surface  velocity  representations  respectively,  choose  a sot  of 
collocation  points,  and  finally  select  an  efficient  numerical  integration 
routine, 

Ifl*  A surface  S*  nay  bo  defined  as  the  continuous  union  of  m sub-regions  S^*, 

each  ol'  which  interpolates  to  some  degree  a set  of  nodes  lying  on  S,  If  the 
transforma ti on 

r = r(u,v)  07) 


maps  a point  (u,v)  in  a local  surface  coordinate  system  for  the  subregion 

to  a position  £ = (x,  y,  s)  in  the  global  cartesian  coordinate  aystom  for  the 
sisrfaoe  S,  then  an  integral  of  the  form 


may  be  approximated  by 

n . 

I w E I f(r)  dudv 

i=1  4 • 

®i 

where 


(38) 


(39) 


(40) 


and 


(41) 


(D^,  Dg,  D^)/v^  is  the  unit  normal  to  at  (u,  v),  irtiloh  is  required  to 
evaluate  the  normal  derivative  of  the  Green's  function. 


42,  In  particular,  triangular  subregions  of  S are  Interpolated  quadra tioally 
through, 

r r (l-u-v)(l-2u-2v)  ♦ u(2u-1 ) ♦ v(2V-1 ) 

(42) 

♦ 4uvr^  ♦ 4»(l-u-v)r^  ♦ 4u(l-u-v)r^ 


wliero  the  nodes  r, , r,  and  r,  map  to  the  vertices  (O,  O),  (l , 0)  and  (0,  1 ) 
respectively  of  a triangle  in  the  (u,  v)  surface  coordinate  system,  and  the 
nodes  ■•P  midside  positions  (^,  ^)  (O,  ^)  and  (i»  o) 

respectively.  To  ensure  an  invertible  napping  with  the  Jacobian  being 


13. 


non-vaniahing  within  the  triangle  certain  reatrlotiona  are  neeeasary  upon  the 

poaitions  of  the  midaide  nodea  and  r^.^^  A aimilar  trana formation  will 

map  biquadratieslly  nine  reaaonabiy  poaitioned  nodea  to  the  vertioea,  midaide 
poaiUona  and  centre  of  a unit  acpare  in  the  local  aurfaoe  coordinate  ayatem. 


43.  In  the  preaent  iB5)lementation  of  the  CHIEP  method  each  aubregion  ia 

thua  defined  by  quadratioally  interpolating  either  aix  or  nine  nodea  lying  in 
S to  give  triangular  or  quadrilateral  aubregiona  reapeotively  aa  appropriate. 
Piecewiae  oonatant  approximationa  are  taken  for  the  aurfaoe  preaaure  and  normal 
aurfaoe  velocity,  auch  that  the  baala  funotiona  are  defined  aa 


; 


r in  S.* 
~ i 

otherwiae 


(43) 


44.  In  the  atandard  finite  element  terminology,  theae  elements  are  of  the 
super parametric  type  with  the  geometry  being  specified  to  a greater  accuracy 
than  the  unknown  functions. 


45.  Hie  collocation  points  are  naturally  taken  as  the  centroids  of  the  surface 
elements  and  the  piecewise  constant  ooii^puted  solution  is  usually  observed  to 
best  ^proximate  the  true  solution  at  these  poaitions.  Althou^  such  a function 
approximation  may  beem  rather  crude,  it  has  successfully  formed  the  basis  of 

many  practical  acoustic  radiation/soattering  computer  programs.^ ^ 

46.  Numerical  experience  with  a spherical  surface  indicates  that  quadratic 
surface  elements  lead  to  an  increase  in  accuracy  over  a piecewise  planar  surface 
by  a factor  of  at  least  20,  whilst  a linear  fVmotion  approximation  gives  little 
or  no  improvement  over  piecewise  constant  basis  funotiona.  The  latter  behaviour 
13  analogous,  for  integral  equations,  to  the  midpoint  rule  being  more  accurate 
than  the  trapeaium  rule  for  numerical  integration.  The  current  interest  in 
Integral  equation  methods  will  no  doubt  provide  more  insight  into  the  gains  to 

be  made  with  higher  order  approximations.  In  particular,  Burton^^  has  suggested 
the  use  of  isoparametric  bicubic  spline  elements  as  the  basis  of  an  improved 
acoustic  program. 

^7.  The  coefficients  of  the  acoustic  matrices  must  be  evaluated  through 
numerical  quadrature.  Singular  integrands  arise  whenever  the  collocation  point 
r^  lies  within  the  region  of  integration  Sj*,  that  is,  for  the  diagonal  matrix 

elements  and  i=1  ,2,  ...,  m derived  from  the  SHB, 


48.  The  dominant  part  of  the  singularity  in  the  diagonal  elements  of  the 
A matrix  may  be  removed  by  rewriting  the  double  layer  potential  term  as 

I = I (p(£)  |~  (r,r‘)  - p(r')  ^ (r,r')|  dS^  (A4) 

-^p(r') 

where 


C^(r,r') 


/I ' ' 


1 

4ff  I i^r ' I 

0^  /w 


(45) 


The  weak 


3G 

aingularlty  of  (£»£*)  bean  aubtraoted  out  by  one  of  the 


same  form  of  t — (r.r')# 
on  ^ ^ 
r 

49*  Then,  for  the  particular  basis  Ainctions  defined  by  (43), 

n bG  p 

*ii  = ’•“  * jf,  i . ^ i % “■■■ 

K* 


30.  Nuaerioally,  the  same  integration  rule  should  be  eiqilpyed  for  both  inte- 
grals (srer  the  subregion  S^*  to  aohiewe  subtaotion  of  the  singularity. 


Unfor^nately,  it  is  not  possible  to  treat  the  weak  sii^ularity  of  the  single 


layer  potential  term  in  a siodlar  manner,  and  a simpler  but  leas  reliable  scheme 
is  employed.  The  region  is  further  subdiTided  into  subregions  suoh  that  the 

collocation  point  lies  at  a vertex  of  each  subregion.  In  the  surface  co- 
ordinate system  (u,  t),  the  collocation  point  elthBr  lies  at  the  positlan  (4, 

in  the  unit  square  or  at  ("j,  -j)  within  the  standard  triangle.  The  unit  square 

^ ^ 11 

is  divided  into  four  subsquares  by  the  lines  u = and  ▼ = '2  standard 

triangle  is  divided  into  three  aubtrlanges  with  the  point  •j)  being  a oomnon 

vertex.  The  integration  rules  adopted  are  Gaussian  rules  with  degree  of  pre- 
cision 3 for  non-singular  integrands  and  degree  of  preolslon  7 within  each  sub- 

region  for  singular  integreuids.^*^^  Highar  order  rules  wore  found  to  give  no 
significant  inox^ase  in  accuracy  whilst  with  lower  ox^er  rules  the  final  solution 
detex^orated  in  accuracy. 


51.  Advantage  may  be  taken  of  any  goons tri cal  symmetry  in  the  body  idien  the 
surface  pressure  and  normal  surface  velocity  have  an  Identical  symmetry.  The 
same  unknown  parameters  may  be  associated  with  a number  of  surface  elements 

and  consequently  it  is  only  necessary  to  oollooate  at  the  centroid  of  one  of 
those  elements.  In  particular  for  complete  axial  symmetxy,  as  in  the  examples 
discussed  below,  this  technique  may  bo  and  is  osployed.  In  suoh  oases  however, 
it  would  be  considerably  more  efficient  to  take  full  advantage  of  the  synmetry 
and  reformulate  the  px*oblem  in  terms  of  a line  rather  than  a surface  integral 

equation.^  ^ 


32.  As  an  Illustration  of  the  numerical  method,  consider  a plane  acoustic  wave 
incident  upon  a submerged  homogeneous  isotropic  spherical  shell  with  no  material 
damping  enclosing  a vaouum.  Two  particular  cases  are  presented;  firstly,  a 
relatively  thick  shell  which  acts  very  much  like  a rigid  soatterer,  and  secondly 
a thin  shell  where  the  stxnicture  has  a more  pronounced  effect  cn  the  acoustic 
field. 


i 
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53. 


Axisymmetric  finite  elements  are  used  to  model  the  shell  for  which  the 

mass  and  stiffness  matrices  have  been  defined  many  times. A six-noded  iso- 
parametric quadrilateral  with  quadratic  variation  in  the  circumferential 
direction  and  linear  variation  in  the  radial  direction  was  found  to  give  accept- 
able results  with  one  element  through  the  thickness  of  the  shell.  The  acoustic 
and  structural  elements  are  illustrated  in  Figure  1,  Structural  nodal  lines  on 
the  outer  spherical  surface  coincide  with  acoustic  element  boundaries  in  order 
to  facilitate  computation  of  the  interaction  matrix  L. 


54,  Firstly,  consider  a spherical  steel  shell  of  outer  raidius  2.0  cm  and 
inner  radius  1.5  cm  submerged  in  water,  Tlie  frequency  of  the  incident  plane 
wave  is  chosen  such  that  the  acoustic  wavenumber  k (=w/c)  is  unity.  The 
scattered  sui'face  pressure  is  initially  computed  assuming  the  outer  surface  of 
the  shell  to  oe  perfectly  rigid,  in  which  case  equations  {^0)  alone  fully  des- 
cribe the  acoustic  field  (with  v a O).  In  Figure  2 the  continuous  lines  are  the 


42 


real  and  imaginary  parts  of  the  exact  solution,  while  the  computed  solution  is 
shown  using  18  elements  in  the  direction  of  0 and  24  elements  in  each  band  around 
the  axis  of  symmetry. 


55.  When  the  response  of  the  elastic  structure  is  taken  into  account,  modelled 
with  18  structural  elements,  the  results  are  as  illustrated  in  Figures  3 and  4* 
The  excict  solution,  shown  as  the  continuous  lines,  is  obtained  tlirough  a separa- 
tion of  variables  analysis  following  Ctoodman  and  Stern. ^ In  Figure  3 the  exact 
solution  for  the  rigid  surface  is  also  shown  for  comparison. 


56.  Now  consider  a thin  spherical  shell  of  outer  radius  5.0  inner  radius 

4.921875  cm.  The  scattered  surfauie  pressure  when  the  outer  surface  is  considered 
to  be  rigid  is  given  in  Figure  5i  using  the  same  number  of  acoustic  elements  as 
for  the  smaller  spliere.  Figures  6 and  7 show  the  effect  of  including  the  vibra- 
tional response  of  the  steel  shell  through  I8  axisymmetric  structural  elements. 
The  exact  rigid  sphere  solution  is  also  superimposed  upon  Figure  6 to  emphasis 
the  stronger  fluid-solid  interaction  in  this  case.  This  is  also  evident  from  a 
comparison  of  the  magnitude  of  the  surface  displacements  in  the  two  cases. 


57.  In  neither  of  these  problems  is  it  necessary  to  consider  interference  from 
interior  eigenfunctions.  For  a spherical  fluid— solid  interface  non— unniueness 
can  only  occur  when  the  non-iiimensional  parameter  ka  (wavenumber  x radius  of 
sphere)  is  a zero  of  the  n-th  order  spherical  Bessel  function  of  the  first  kind. 


for  some  value  of  n. 


CONCLUSIONS 
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58.  By  combining  equations  derived  from  a finite  element  analysis  of  a vibrating 
elastic  structure  with  tliose  from  an  integral  equation  representation  of  lui 
infinite  exterior  acoustic  field,  an  analysis  of  the  coupled  dynamic  interaction 
problem  has  been  shown  to  be  fcat'ible. 


59.  The  structural  and  ai'oustic  models  may  be  chosen  completely  independently 
although  the  procedure  does  simplify  if  interpolation  luxies  are  made  to  coincide, 
the  coupling  being  defined  through  three  interaction  matrices  X,  L and  L'.  Thus 
advanteige  may  be  taken  of  existing  computer  program  packages  designed  to  solve 
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the  two  uncoupled  problems.  In  particular,  it  is  hoped  to  combine  the  PAPEC 
structural  analysis  program  with  the  acoustic  radiation/scattering  program  des- 

cribed  in  the  previous  section. 


60,  borne  test  examples  for  which  analytic  solutions  are  available  have  demon- 
strated the  practicability  and  accuracy  of  the  ^proach.  It  is  hoped  in  the  near 
future  to  be  able  to  compare  numerical  euid  experimental  results  for  some  more 
complex  structures  in  water.  If  an  unacceptably  large  numer  of  degrees  of  free- 
dom for  the  structural  model  are  reqpiired,  it  may  then  be  necessary  to  consider 
the  alternative  modal  approach. 
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